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Abstract. In this paper we focus on the TP-AGB evo- 
lution of intermediate-mass stars experiencing envelope 
burning (M = 4-=-5M ). Our model of the TP-AGB phase 
is suitably designed to follow the peculiar behaviour of 
these stars, to which the simple analytical treatment valid 
in the low-mass range can no longer be applied. 

The approach we have adopted is a semi-analytical one 
as it combines analytical relationships derived from com- 
plete models of TP-AGB stars with sole envelope models 
in which the physical structure is calculated from the pho- 
tosphere down to the core. The solution for the envelope 
models stands on an original numerical method which al- 
lows to treat major aspects of envelope burning. 

The method secures that, during the quiescent inter- 
pulse periods, fundamental quantities such as the effective 
temperature, the surface luminosity, the physical struc- 
| ture of the deepest and hottest layers of the envelope, 
and the related energy generation from nuclear burning, 
are not input parameters but the consequence of envelope 
model calculations. This minimizes the use of analytical 
relations, thus giving our results greater homogeneity and 
accuracy. 

Moreover, we would like to draw the attention on the 
general validity of our algorithm which can be applied also 
to the case of low-mass stars, in which envelope burning 
does not occur. 

Our efforts are directed to analyse the effects produced 
by envelope burning, such as: i) the energy contribution 
which may drive significant deviations from the standard 
core mass-luminosity relationship; and ii) the changes in 
the surface chemical composition due to nuclear burning 
via the CNO cycle. 

Evolutionary models for stars with initial mass of 
4.0, 4.5, 5.0-Mq and two choices of the initial chemical corn- 
Send offprint requests to: P. Marigo 



position ([Y = 0.28, Z = 0.02] and [Y = 0.25, Z = 0.008]) 
are calculated from the first thermal pulse till the complete 
ejection of the envelope. We find that massive TP-AGB 
stars can rapidly reach high luminosities (—6 > Mb ] > 
—7), without exceeding, however, the classical limit to the 
AGB luminosity of Mboi — —7.1 corresponding to the 
Chandrasekhar value of the core mass. No carbon stars 
brighter than Mboi ~ —6.5 are predicted to form (the 
alternative of a possible transition from M-star to C-star 
during the final pulses is also explored) , in agreement with 
observations which indicate that most of the very lumi- 
nous AGB stars are oxygen-rich. 

Finally, new chemical yields from stars in the mass 
range 4 5M are presented, so as to extend the sets of 
stellar yields from low-mass stars already calculated by 
Marigo et al. (1996). For each CNO element we give both 
the secondary and the primary components. 

Key words: stars: evolution - stars: AGB and post- AGB 
- nuclear reactions, nucleosynthesis, abundances - stars: 
abundances - stars: carbon - ISM: Planetary Nebulae: 
general 



1. Introduction 

This study deals with the theoretical evolution of TP- 
AGB stars undergoing envelope burning - i.e. those in 
which the base of the convective envelope deepens into 
high temperature regions so that H-burning via the CNO 
cycle can occur. This circumstance is usually encoun- 
tered in AGB stars with massive and extended envelopes 
(M cnv 2b 2.5M©), as indicated by several authors (Sug- 
imoto 1971; Uus 1972; Iben 1973; Renzini & Voli 1981; 
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Scalo et al. 1975; Blocker & Schonberner 1991; Lattanzio 
1992; Boothroyd & Sackmann 1992; Vassiliadis & Wood 
1993; Blocker 1995a; Marigo et al. 1996). The most no- 
table signature of the TP-AGB evolution of these stars 
is the break-down of the core mass-luminosity relation at 
high luminosities (Blocker & Schonberner 1991; Lattanzio 
1992; Boothroyd & Sackmann 1992). 

Previous analyses of the AGB phase performed with 
the aid of envelope models have already pointed out that 
nuclear burning can take place in the deepest and hottest 
layers of the convective envelope (Scalo et al. 1975; Renzini 
& Voli 1981; Marigo et al. 1996). 

In Marigo et al. (1996) we adopted the assumption 
that those regions do not provide energy to the star, and 
the constraint that the core mass-luminosity relation is 
always satisfied. 

The notable improvement of the present approach is to 
account for the energy generated by the burning regions of 
the envelope. For this purpose, we develop an algorithm 
to include nuclear burning in the envelope of a star in 
thermal equilibrium. 

Establishing how deep in the interior the convective en- 
velope penetrates is crucial to determine the efficiency of 
the nucleosynthesis processes, for the extreme sensitivity 
of nuclear rates to temperature and density. Furthermore, 
nuclear burning changes not only the chemical composi- 
tion of the envelope, but also affects the luminosity, the ef- 
fective temperature, and the radius of the star. Therefore, 
the calculation of the envelope models must be carried out 
with great accuracy. 

The theoretical approach developed in this study of- 
fers some advantages: it keeps the agility of purely an- 
alytical models for quick computing and easy testing of 
different prescriptions but it has greater accuracy and self- 
consistency, which draw it closer to full calculations of the 
TP-AGB phase. 

This paper is organized in three main sections. In the 
first part (Sect. 2) we summarize the physical conditions 
under which the core-mass luminosity relation is found 
and we analyze the break-down of this due to envelope 
burning contributing to the total stellar luminosity over 
the quiescent inter-pulse period (Sect. 3). The theoretical 
method developed to this aim is described and discussed 
(Sect. 3.1). 

In the second part (Sect. 4) we apply our method to 
follow the TP-AGB phase of the 4.0, 4.5, 5.0M Q stars with 
two choices of the initial chemical composition, namely 
[Y = 0.25, Z = 0.008] and [Y = 0.28, Z = 0.02]. The 
calculations are carried out from the end of the E-AGB 
up to the complete ejection of the envelope. 

In the third part (Sect. 5) we study the nucleosynthesis 
associated to envelope burning which, together with the 
third dredge-up, concurs to alter the surface abundances 
of the elemental species. Finally (Sect. 6), we provide new 
chemical yields from stars in the above mass range, distin- 



guishing for each element of the CNO group the primary 
and secondary contribution. 

2. The core mass-luminosity relation 

Paczynski (1970) first discovered the existence of a sim- 
ple, almost linear relationship between the core mass of 
a double-shell burning star and its quiescent inter-flash 
luminosity along the TP-AGB. 

Over the years, continuous upgrading of the input 
physics (e.g. new opacities, revised nuclear reaction rates, 
equations of state, etc.) and detailed calculations of the 
AGB phase for wide ranges of stellar masses and chemical 
compositions resulted into a flourishing of different core 
mass-luminosity (M c — L) relationships (Ibcn 1977; Ibcn 
& Truran 1978; Havazelet & Barkat 1979; Wood & Zarro 
1981; Lattanzio 1986; Boothroyd & Sackmann 1988a). 

In this work (see also Marigo et al. 1996) we adopt the 
following relations: 

L Mc = 238 000^ 3 ZgSo(M c 2 - 0.0305M C - 0.1802) (1) 

for stars with core mass in the range 0.5M Q < M c < 
O.66M (Boothroyd & Sackmann 1988a); and 

L Mc = 122585Ai 2 (M c - 0.46)M al9 (2) 

for stars with core mass M c > O.95M (originally from 
Ibcn & Truran 1978 according to the slightly modified 
version of Groenewegen & de Jong 1993; see also Ibcn 
1977). 

In these formulas stellar masses M and luminosities L 
are in solar units; Zqno is the total abundance (in mass 
fraction) of C, N, O isotopes in the envelope; \x = 4/ (5A + 
3 — Z) is the mean molecular weight under the assumption 
of a fully ionized gas, where X and Z are the abundances 
of hydrogen and metals, respectively. For stars with core 
mass 0.66M Q < M c < 0.95M Q a linear interpolation is 
performed. 

In our model we also account for the steep luminosity 
increment as a function of the core mass during the first 
inter-pulse periods, before the onset of the "full ampli- 
tude" regime to which the standard M c — L relation given 
by Eqs. (1) or (2) refers. The initial subluminous evolution 
is followed by adopting a simple linear relation between 
the luminosity and the core mass, Lfi rs t = S x M c + Q, 
with the slope S given as a function of the initial mass M; 
of the star: 

S = 60761.2 x exp(M;/2.) (3) 

where masses and luminosities are expressed in solar units. 
This relation is a least-square fit to the results presented 
by Vassiliadis & Wood (1993) (see their Fig. 12). For each 
value of Mi, the intercept Q is fixed by the values of the 
luminosity and core mass at the first thermal pulse. 

The existence of the core mass-luminosity relation for 
AGB stars was given a transparent analytical explanation 
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by Tuchman et al. (1983) just starting from the equations 
of stellar structure under specific physical conditions ( "ra- 
diative zero solution" ; sec Egglcton f 967; Paczynski f 970). 
The stellar structure is required to possess a degenerate 
core of mass M c surrounded by a narrow radiative shell 
(or double shell) source providing the luminosity beyond 
which there must exist a thin (with a mass AM <§C M c ) 
and inert (the luminosity is constant) transition region in 
radiative equilibrium. This is bounded at the top by the 
base of the convective envelope. Integrating the equations 
of hydrostatic and radiative equilibrium over the tran- 
sition region, after simple substitutions, Tuchman ct al. 
(1983) derived an analytical expression relating the lu- 
minosity L of the star to its core mass M c , with some 
dependence on the chemical composition and the temper- 
ature just above the burning shell. The effect of the outer 
envelope on this relation is negligible because of the ex- 
tremely steepness of pressure, temperature, and density 
gradients in the radiative zone. Thus, during the quies- 
cent TP-AGB phase, the evolution of the core is expected 
to be essentially decoupled from that of the envelope (no 
mass dependence). 

However, Tuchman ct al. (1983) suggested that the 
M c — L relation is likely to break down for stars with 
envelope burning, as the base of the convective envelope 
may penetrate so deeply as to eat up the radiative buffer 
over the H-burning shell. 

Such a possibility was anticipated by Scalo et al. (1975) 
who pointed out that, during the quiescent inter-pulse pe- 
riod, in massive and luminous AGB stars the envelope 
convection may extend all the way down to the H-burning 
shell, as a result of the high radiation pressure Pr. In brief, 
when the contribution of Pr to the total pressure becomes 
significant, the radiative (V ra d) and adiabatic (V a d) tem- 
perature gradients are expected to decrease less steeply 
with depth, and to reach the minimum value of 0.25, re- 
spectively. Both facts concur in drawing the base of the 
convective envelope - reached when V ra d = V a d - deeper 
in the interior. High temperatures (7b > 40 x 10 6 K) can 
be attained, thus allowing nuclear burning at the base of 
the envelope. 

3. Envelope burning as an energy source 

Over the years it has become clear that the energy genera- 
tion due to envelope burning can make massive TP-AGB 
stars to significantly deviate from the M c — L relation. 
High luminosities are soon reached (Mboi < —6), at least 
until the residual envelope is massive enough to allow for 
high temperatures (T b > 40 x 10 6 K) at its bottom. All 
this is indicated by complete models of TP-AGB stars 
(Blocker & Schonberner 1991; Lattanzio 1992; Boothroyd 
& Sackmann 1992; Vassiliadis & Wood 1993). 

Earlier analyses, primarily aimed to study the struc- 
ture of envelope burning and associated nucleosynthesis, 
were carried out with the aid of convective envelopes mod- 



els (Scalo et al. 1975; Renzini & Voli 1981). The same kind 
of theoretical approach adopted in the semi-analytical TP- 
AGB model by Marigo et al. (1996). The assumption com- 
mon to all those models is that the envelope satisfies the 
classical core mass-luminosity relation. This is no longer 
valid in massive AGB stars, due to the energy genera- 
tion by envelope burning. Moreover, numerical integra- 
tions were made under the constraint that the radial co- 
ordinate of the base of the convective envelope coincides 
with the core radius (Renzini & Voli 1981; Marigo et al. 
1996). This approximation can be risky in the case of mas- 
sive TP-AGB stars as the base of the convective envelope 
falls in a region, which is very thin in mass but across 
which the structural gradients are very steep. Therefore, 
fixing a priori the inner extension of the envelope convec- 
tion may lead to heavily under- or over-estimate the base 
temperature, to the value of which the nuclear rates are 
highly sensitive. In contrast, the correct treatment of the 
envelope burning requires that the energy being generated 
by nuclear reactions is taken into account and the exact lo- 
cation of the deepest layer through which external convec- 
tion penetrates is determined (e.g. by the Schwarzschild 
condition). 

3.1. The method of solution 

First of all, we relax the original assumption fixing the 
local luminosity at a constant value (equal to L) through- 
out the envelope (i.e. no energy source). This condition is 
replaced with the equation of energy balance: 

^ = 47rr 2 e r p r (4) 
dr 

where L r is the local luminosity; p v the gas density; e r is 
the amount of energy released by nuclear transmutation 
of hydrogen into helium per unit of mass and time. All the 
variables are evaluated at the radius r. 

The integration of Eq. (4) across the envelope must 
satisfy the basic constraint given by the principle of en- 
ergy conservation. This means that during the quiescent 
inter-pulse period the surface luminosity L is the sum of 
different contributions, namely those from: the radiative 
H-burning shell, Lr; the nuclear burning at the base of the 
convective envelope, Leb; the small contributions from the 
He-burning shell, Ln e , the gravitational contraction of the 
core, Lq, and the rate of energy loss via neutrinos, L v : 

f Rs dL 

L= —-^dr = L G + Lu c + Lu + L BB -L u (5) 
Jo dr 

where Rs = \fL] (47rerT c 4 ff ) is the stellar surface radius. In 
this study we do not consider neutrino losses. This nega- 
tive contribution is very small compared to the total lumi- 
nosity, and it is indeed negligible in low-mass stars (Ibcn 
& Tutokov 1984; Koster & Schonberner 1986). In massive 
cores (i.e. M c ^ 0.8M©), its absolute value may be ini- 
tially comparable to the gravitational energy release from 
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core contraction during the quiescent inter-pulse period 
(Blocker 1995a, b). However, as the star evolves the core 
interiors become cooler and more degenerate so that the 
neutrino luminosity falls again well below the gravother- 
mal luminosity. 

Each term (but for L v ) of the sum in Eq. (5) can be 
explicitly expressed as: 



Lg + Luc 
L H 





dL T 
——dr 


Jo 


ar 
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——ar 




dr 


rRs 


dL T 
——dr 


J Rconv 


dr 



(6) 



where i? coro is the core radius, here defined as the radius 
at the bottom boundary of the H-burning shell, and R conv 
is the radial coordinate of the base of the convective en- 
velope. In the system of equations (6) the only quantity 
that can be directly evaluated with our envelope model is 
Leb, whereas for the other terms we need some analytical 
prescriptions. We make use, in particular, of the relation 



L G + L Hc = 2000(M/7) U 4 exp[3.45(M c - 0.96)] 



(7) 



derived from the results of complete calculations of the 
AGB phase (Iben 1977; see also Groenewegen & de Jong 
1993). 

As a first approximation, we consider the energy con- 
tribution from envelope burning, Leb, as a correction term 
which is added to the luminosity Lm c predicted by the 
M c — L relation, Eqs. (1) and (2), to get the total lumi- 
nosity, i.e. L ~ Lm c + Leb- Substituting L into Eq. (5), 
we can express the luminosity provided by the radiative 
H-burning shell as: 



L„ = L 



M, 



L G 



(8) 



There is one aspect of the above scheme deserving some 
remarks. We have assumed that Lh and Leb are decou- 
pled, whereas in a real star with envelope burning as soon 
as the envelope convection penetrates into the H-burning 
region the two luminosities are related to each other in 
a way that cannot be ascertained a priori without cal- 
culating a complete stellar model. A possible way out of 
this difficulty could be a scheme in which the structure 
of the radiative double-shell source with "deep envelope" 
is integrated down to the carbon-oxygen core. This ap- 
proach could be explored in a future study. On the other 
hand, a preliminary analysis shows that the relative error 
of our approximation for Lh (see Eq. (8)) with respect to 
the results from full AGB calculations is mostly comprised 
within ~ 2% 10% (Marigo 1998, in preparation). 

However, despite the adopted approximation, the 
present model is a step forward when compared with pre- 
vious analyses of envelope burning performed with the 
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Fig. 1. Solution of the system of two equations (12) in the 
plane L — T e g. The solid line corresponds to the locus where 
F(L,T e ff) = (L(R conv ) — Lm c ); the asterisk marks the point 
(solution) where G(L,T eff ) = (Q C orc = Q(l Hc +l g ))- 



aid of purely analytical assumptions (Groenewegen & de 
Jong 1993), or static envelope models (Marigo et al. 1996). 
Moreover, the results obtained with this scheme turn out 
to be consistent with those from full evolutionary calcula- 
tions. 

3.2. Boundary conditions 

Given the stellar mass M, the core mass M c , the chemi- 
cal composition of the envelope, the luminosities Lh, Lne, 
Lg at a certain time, the task is to find T e g and Leb 
(and hence total luminosity L) of the stellar model sat- 
isfying Eq. (5). The solution is univocally singled out by 
means of envelope integrations constrained by two bound- 
ary conditions, which naturally derive from the system of 
equations (6) : 



L(R C 
L(i? c( 



= L G 
= L G 



Lhc 
LHe 



Li 



(9) 
(10) 



We define the variable Q(r) = M(r)/M as the mass 
coordinate of the layer of radius r, where M(r) is the 
mass contained below it and M is the total stellar mass. 
Let us call Q corc and Q(l Ho +l g ) tne mass coordinates of 
the core and of the layer with a luminosity Lh + Lq, 
respectively. According to the adopted definition of the 
core, it is physically equivalent to replace Eq. (9) with: 



<9(L H c + Lg) 



(11) 
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Fig. 2. Physical conditions in the deepest layers of the model 
corresponding to the case illustrated in Fig. 1. The variables 
V r ad, Vad, P = fgas / P , L r /L are plotted as function of 
the radial coordinate. The arrows indicate the width of the 
H-burning shell, with the upper boundary defined by an en- 
ergy generation rate of 10 erg gr _1 s _1 . A few values of the 
fractional luminosity L T /L are displayed at interesting points 
(marked by full circles), namely: i) the inner boundary of the 
H-burning shell at -R C or e ; ii) the base of the convective enve- 
lope (when Vrad = Vad, asterisk) at R con v', and iii) the level 
at which the local luminosity amounts to the 99% of the total 
luminosity. In this case the envelope convection extends deeply 
in the interior of the H-burning shell. 

which is more suited to our calculations. To summarize, we 
deal with two boundary conditions (Eqs. (10) and (11)), 
which require to solve a system of two equations in the 
two unknowns L and T e g: 

(G(L,T cS ) = Q COIC - Q(l Ho+ l g ) = . . 

\F(L, T off ) - L(R conv ) - (L G + L Hc + L k ) = [LZ > 

The adopted numerical technique is Brent's method 
for root finding (Numerical Recipes 1990) applied to both 
functions. A tolerance er, = 10~ 5 on the root of function F 
(expressed as the difference of logarithmic luminosities), 
seems reasonable for our purposes. Because of the extreme 
steepness of the structural variables across the H-burning 
shell, we need to work with a high accuracy, eq = 10~ 6 , on 
the root of function G, such that iterations will continue 

Until |Qcoro - Q(L H c + Lg)| < e Q- 

This method proves to be general and able to work also 
when envelope burning is not operative, e.g. in the final 
stages of high mass stars when the mass of the residual 
envelope is drastically reduced, or in the whole TP-AGB 
evolution of low-mass stars which never deviate from the 



Fig. 3. The same as in Fig. 2, but referred to the stage soon 
after the end of the E-AGB phase. In this case a radiative 
transition region exists between the upper boundary of the 
H-burning shell and the base of the convective envelope so 
that the core mass-luminosity relation still holds. 

M c — L relation. In these cases Leb = so that the local 
luminosity all over the envelope coincides with the surface 
luminosity given by L = Lm c ■ The only unknown quantity 
is then T e g which is the root of the function G(Lm c , T c g). 

In order to show how the method works, we con- 
sider the case of a stellar model with current mass M = 
4.9226M©, core mass M c = 0.9502M Q , and luminosity 
logLM c = 4.5175, this latter being derived from the M c —L 
relation. These values refer the final stage of the 20 th 
inter-pulse period of a 5M Q , [Y = 0.25, Z = 0.008] star. 
In the log L — log T c g plane of Fig. 1 we show the locus 
(solid line) of the models satisfying the condition F = 
of Eq. (10). Each i th point along the line is derived from 
envelope integrations carried out with an input luminosity 
Li = (Lq + Lfjo + ^h) + ALj, where ALi is an arbitrary 
amount of extra-luminosity. At given Li there is a unique 
value of the effective temperature, T e g y i (hence a unique 
envelope model) root of the function F. Along the right- 
hand side vertical axis, several values of the function G 
- e.g. G(Li,T e ff^i) - derived from envelope models satis- 
fying the condition F = are displayed in a one-to-one 
correspondence with the luminosity. 

The solution of the system of Eqs. (12) is represented 
by the pair of values (logL = 4.62342, logT cff = 3.52689) 
at which both F = and G = 0. This is indicated by 
the asterisk along the solid line. In this particular model, 
the contribution from envelope burning, Leb, amounts to 
~ 22% of the total luminosity. 
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Figure 2 shows a few physical quantities characterizing 
the structure of this model in the deepest regions of the 
envelope and beyond, down to the lower boundary of the 
H-burning shell, i? C oro- To be noticed is that the base of the 
convective envelope is reached in the deep interior of the 
H-burning shell, where the temperature is lb ~ 100 x 10 6 . 

For comparison, the same quantities are plotted in Fig. 
3, but relative to the 5M Q , [Y = 0.25, Z = 0.008] star soon 
after the end of the E-AGB phase, with M c = 0.9344M o , 
logT cff = 3.5583, and L Mc = 4.3502L©. At this stage, 
envelope burning is very weak (Tb ~ 24 x 10 6 ; Leb = 0), 
so that the star still obeys the M c — L relation. Note the 
existence of the radiative buffer between the H-burning 
shell and the base of the convective envelope. 



4. Evolutionary calculations 

With the method described in Sect. 3.1 we calculated 
the TP-AGB evolution of stars with initial mass of 
4.0, 4.5, 5.0M Q and two choices of the chemical compo- 
sition [Y = 0.28, Z = 0.02] and [Y = 0.25, Z = 0.008]. 
In the following we briefly recall the basic ingredients 
adopted in the model, the same ones employed by Marigo 
et al. (1996) to whom the redear should refer for more 
details. 

The initial conditions at the first thermal pulse are 
taken from full evolutionary calculations extending from 
the main sequence till the beginning of the TP-AGB phase 
kindly provided us by Leo Girardi (1997, private commu- 
nication) . For each stellar mass we single out the first sig- 
nificant pulse, usually defined as the first He-shell flash 
in which L^ x > L, where L is the surface luminosity 
(Boothroyd & Sackmann 1988b). All the quantities of in- 
terest are then referred to the time corresponding to the 
quiescent pre-flash luminosity maximum in the light curve. 
This procedure let us account for the changes in the sur- 
face chemical composition caused by the previous mix- 
ing processes, namely, the first dredge-up as the star first 
reaches its Hayashi track after central hydrogen exhaus- 
tion, and the second dredge-up at the base of the E-AGB 
(Tables 2 and 3). 

It is worth noticing that the input physics of our enve- 
lope model is the same as in the Padua stellar evolution 
code, so that our results on the TP-AGB phase are homo- 
geneous with the previous evolutionary history. The main 
characteristics of these stellar models can be summarized 
as follows: 

- The stellar models in usage here to derive the ini- 
tial conditions allow for convective overshoot from the 
core and external envelope (Alongi et al. 1993; see also 
Marigo et al. 1996). Therefore, the maximum initial 
mass, M up , for a star to build an electron-degenerate 
C-0 core, i.e. to evolve through the AGB phase, is 
around 5M Q . Classical models without overshoot pre- 



— The high-temperature opacities are those of Iglesias 
& Rogers (1996). They are particularly suitable to 
AGB calculations inclusive of the third dredge-up as 
the authors also provide tables for a variety of carbon- 
and oxygen-rich mixtures. At low temperatures (T < 
10 4 K) we use the molecular opacities of Alexander & 
Ferguson (1994). 

— In the outermost super-adiabatic convective region, 
the adopted value of the mixing length parameter is 
a = 1.68, according to the calibration of the solar 
model. 

— Finally, the CNO-cycle nuclear reaction rates are from 
Caughlan & Fowler (1988). 

The analytical relations used to construct our TP- 
AGB synthetic model are: 

(1) The core mass-luminosity relation, given by Eq. (1) 
or Eq. (2). 

(2) The core mass-interpulse period relation, taken 
from Boothroyd & Sackmann (1988b): 



logT;, 



4.50 (1.689 -M c ) for Z = 0.02 
4.95 (1.644 -M c ) for Z = 0.001 



(13) 



where Tj p and M c are expressed in years and solar units, 
respectively. 

(3) The rate of evolution, according to Groenewegen 
& de Jong (1993): 



^ = 9.555 10- 12 ^ 
dt X 



(14) 



where dM c /dt is the growth rate of the core mass 
(Moj/r -1 ), X is the hydrogen abundance (in mass frac- 
tion) in the envelope, Lh is the luminosity produced by 
H-burning (in solar units). The numerical factor takes into 
account the energy released from the nuclear conversion 
of 1 g of hydrogen into helium (~ 6.4 10 18 erg). 

(4) Flash-driven luminosity variations - We apply 
proper corrections as a function of the envelope mass 
(Groenewegen & dc Jong 1993) to account for the de- 
viations in the luminosity caused by the occurrence of 
thermal pulses, i.e. the post-flash luminosity peak and the 
low-luminosity dip. 

(5) Mass loss - The adopted prescription for the mass- 
loss rate is the semi-empirical one of Vassiliadis & Wood 
(1993), which gives the rate as a function of the pulsational 
period of variable AGB stars: 



logM 

M = 6.07023 10~ 3 



11.4 + 0.0123P 
L 



cv, 



(15) 
(16) 



cxp 



diet a higher value, i.e. M up ~ 8M . 



Here, M is given in units of M yr _1 , the stellar luminos- 
ity L is expressed in L , the pulsation period P in days, c 
is the light speed (in km s -1 ) and w OX p (in km s _1 ) denotes 
the terminal velocity of stellar wind. During calculations, 
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at any time step, the mass-loss rate is chosen to be the 
minimum value between those given by Eq. (15) and (16). 
The wind expansion velocity v eX p (in km s^" 1 ) is estimated 
as a function of the period of pulsation: 

v cxp = -13.5 + 0.056P (17) 

with the additional constraint that w cxp lies in the range 
3.0 — 15.0 km s -1 , the upper limit being the typical termi- 
nal velocity detected in high mass- loss rate OH/IR stars. 

The pulsation period P is derived from the period- 
mass-radius relation (Eq. (4) in Vassiliadis & Wood 1993), 
with the assumption that variable AGB stars are pulsating 
in the fundamental mode: 

log P = -2.07+ 1.94 logP-0.9 logM (18) 

where the period P is given in days; the stellar radius R 
and mass M are expressed in solar units. 

For each stellar mass the calculations are carried out 
starting from the first significant pulse until the envelope 
mass, M onv reduces to about O.OO1M that hydrodynam- 
ical studies suggest to be the critical value at which the 
pulsational instability (and hence strong mass-loss) ceases 
(Tuchman ct al. 1979). It is worth recalling that full evo- 
lutionary calculations show that a model star moves off 
the AGB when the residual envelope mass is diminished 
below O.OIMq (Schonberner 1983). From the present re- 
sults it turns out that the reddest point in the H-R dia- 
gram (see Fig. 8) is reached when M env ~ 0.8 ~ O.9M 
in all cases studied, with a typical value logT c ff ~ 3.48 
for models with Z = 0.008, and logT off ~ 3.44 for mod- 
els with Z = 0.02. Then, the tracks slowly move toward 
the blue part of the H-R diagram, lying in the interval 
[3.44,3.48 + 0.03] until M env ~ 0.3 + 0.4M©. The final 
evolution of the models is characterized by a quick de- 
crease of T e g (because of the very small envelope mass). 
However, we do not pretend here to handle the delicate 
topic about the transition from the AGB to the forma- 
tion of planetary nebulae, which is beyond the aim of the 
present study. 

The total number of thermal pulses included in each 
evolutionary calculation is iVp ot = 174 for the [4.0M Q , Z = 
0.008] model; 7V P ot = 146 for the [4.5M ,Z = 0.008] 
model; 7V P ot = 129 for the [5.OM ,Z = 0.008] model. A 
smaller number of thermal pulses are expected to be suf- 
fered by stars with the same initial mass, but lower (e.g. 
solar) mctallicity. Our calculations yield Np ot = 39 for the 
[4.0M Q , Z = 0.02] model; iV P ot = 56 for the [4.5M , Z = 
0.02] model; iV P ot = 83 for the [5.OM ,Z = 0.02] model. 

4-1. Evolution in the M c — L plane 

As already mentioned, the energy contribution from en- 
velope burning makes massive TP- AGB stars to leave the 
M c — L relation during part of their evolution. This effect 
is illustrated in Figs. 4 and 5. During the initial stages of 



the thermally pulsing regime the luminosity evolution is 
governed by the H-burning shell, just re-ignited after the 
second dredgc-up. As the efficiency of the envelope burn- 
ing increases, the star can quickly reach higher and higher 
luminosities (steeply rising parts of the curves), with a rate 
much greater than expected from the slope of the M c — L 
relation (dashed lines). The luminosity growth goes on as 
long as the envelope remains massive enough to allow high 
temperatures at the base of the envelope (see Sect. 4.2 
below). When mass loss becomes significant (and enve- 
lope burning weakens) the luminosity starts to decline and 
the star eventually approaches again the M c — L relation, 
where it lies till the end of the AGB phase. This kind of 
luminosity evolution agrees with the results from full evo- 
lutionary calculations of the TP- AGB phase (Boothroyd 
& Sackmann 1992; Vassiliadis & Wood 1993). 

The above trend may of course vary depending on sev- 
eral factors, part of which are intrinsic to the stars, e.g. 
initial mass and chemical composition, and part to the 
theoretical model, e.g. mixing formalism, mass-loss pre- 
scription, opacities, and nuclear reaction rates (see Rcn- 
zini & Voli 1981; Sackmann & Boothroyd 1991; Boothroyd 
& Sackmann 1992, 1993; Blocker 1995a). 

Figures 4 and 5 show the effect on the above trend 
of mass and chemical composition. The deviation from 
the luminosity predicted by the M c — L relation is more 
pronounced in more massive stars at fixed metallicity. 
Let us consider the [Y = 0.25, Z = 0.008] case (Fig. 4). 
The maximum relative contribution of envelope burning 
to the total luminosity Lbb/L goes from about 0.02 for 
the 4.OM model, 0.14 for the 4.5M Q model, up to 0.23 
for the 5M model. Concerning this latter, the plot shows 
how the decline in luminosity depends on the current value 
of the stellar mass (indicated along the curve), as it is 
progressively reduced by stellar winds. At higher mctal- 
licities (Fig. 5), the overluminosity generated by the enve- 
lope burning is much smaller at given initial mass. In the 
[Y = 0.28, Z = 0.02] case, the maximum relative contri- 
bution is negligible for the 4.OM model (which actually 
never deviates from the M c — L relation), and very small 
for the 5.0 and 4.5M stars, where Leb/L amounts to 
0.15 and 0.05, respectively. 

Concluding this section, we like to briefly comment on 
the dependence of M c — L relation given by Eq. (2) on the 
current star mass (i.e. the factor M 019 ). Ibcn (1977) orig- 
inally introduced a factor M 0A to fit the results for the 
AGB evolution (at constant mass) of his 7.OM star. Sub- 
sequently Iben & Truran (1978) reduced the exponent to 
0.19 to obtain a compromise with the results of Paczynski 
(1970). 

This dependence is indeed quite weak, and it has in 
practice no effect for most of the TP- AGB lifetime. How- 
ever, after the onset of the super-wind, envelope strip- 
ping gets so efficient that the weight of the mass factor 
in Eq. (2) becomes important. In Fig. 4 the final, slow de- 
crease in luminosity at increasing core mass during the last 
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Fig. 4. Luminosity evolution of 4.0, 4.5, 5.0 M© TP-AGB stars 
with composition [Y — 0.25, Z = 0.008] as a function of the 
core mass. The dashed lines correspond to the standard core 
mass-luminosity relation, the filled circles refer to the actual 
stellar luminosity, at its maximum value before each pulse. The 
overluminosity produced by envelope burning rapidly reaches 
a maximum value, it starts to decrease after the onset of the 
super-wind (the numbers along the 5M curve indicate the 
current stellar mass in solar units), finally it vanishes and the 
star approaches again the core mass-luminosity relation. 



z=o.oa 



4xl0 4 




10 4 



0.85 



0.9 0.95 
M core /M e 



Fig. 5. The same as in Fig. 4, but with [Y = 0.28, Z = 0.02]. 



evolutionary stages is just caused by the drastic reduction 
of the stellar mass. 

This trend seems to contradict the existence of the 
core mass-luminosity relation itself, because this latter im- 
plies that the luminosity should not depend on the enve- 
lope mass. We may argue that the factor M 019 actually 
masked some low-efficiency envelope burning (the adopted 
mixing length parameter was a = 0.7 in Iben's calcula- 
tions, roughly corresponding to a ~ 1.4 of other groups). 
Though a cautionary remark on the use of such a relation 
is worthy since it may depend on model details, the gen- 
eral applicability of the our solution scheme is not affected, 
and the adoption of different M c — L relations could be 
explored in a future study. 

4-2. The base temperature 

The temperature at the base of the convective envelope, 
lb, is the key-quantity driving the strength of envelope 
burning. Figures 6 and 7 show T> as a function of the prc- 
fiash luminosity maximum for stars with [Y = 0.28, Z = 
0.02] and [Y = 0.25, Z = 0.008]. Both diagrams share a 
common trend, regardless of the initial mass and chemical 
composition (metallicity) . Starting from the first pulse, as 
the star evolves to higher luminosities a rapid increase in 
Tb occurs in all models. This feature reflects the progres- 
sive inward penetration of the base of the convective enve- 
lope, which is essentially driven by the effect of radiation 
pressure on V ra d and V a( j (see Sect. 2). It is worth notic- 
ing that in these initial stages, when envelope burning is 
growing but still weak, there is a well-defined almost lin- 
ear relation between logTb and logL, with a slope nearly 
independent of the stellar mass. 

The subsequent bending of the Tb curves happens at 
a certain "threshold luminosity" when the envelope con- 
vection eats up the radiative buffer and extends into the 
H-burning shell (Scalo et al. 1975). This critical luminos- 
ity marks the point at which the star begins to deviate 
from the M c — L relation (Figs. 5 and 4). At even higher 
luminosities, the modest increase in Tb is controlled by 
the steep gradient in L T (and V ra d in turn) toward the in- 
nermost layers of the burning region. Once the luminosity 
has reached a maximum value and begins to decline, Tb 
lowers as well because of the reduction of envelope mass 
by stellar winds, and finally it drops so dramatically that 
envelope burning extinguishes. 

Comparing the curves in Figs. 6 and 7 for a given ini- 
tial mass it turns out that, at the same luminosity, a lower 
metallicity favours higher values of Tb along the TP-AGB. 
A similar effect takes place on the characteristic temper- 
ature at which tracks of different metallicities begin to 
flatten out, with a typical value logTb ~ 7.73 for the 
[Y = 0.28, Z = 0.02] case and logT b ~ 7.79 for the [Y = 
0.25, Z = 0.008] case. Moreover, it is worth remarking that 
the lower the stellar mass, the lower is the value of the 
"threshold luminosity". The 4.0M Q , [Y = 028. Z = 0.02] 
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Fig. 6. Temperature T b at the base of the convective en- 
velope as a function of the pre-flash luminosity maximum 
(filled circles) for intermediate-mass TP-AGB stars with 
[Y = 0.28, Z — 0.02]. 

model is a transition case for the occurrence of envelope 
burning. The flattening of the T, curve does not show up 
because as soon as the star reaches the "threshold lumi- 
nosity" (logTb ~ 7.73), the residual envelope is not mas- 
sive enough to sustain high temperatures at its base. As 
a consequence of it, the Tb curve abruptly falls down. 

4-3. The H-R diagram 

Figure 8 shows the evolutionary tracks in the H-R diagram 
for our TP-AGB stars with composition [Y = 0.25, Z = 
0.008] (solid lines) and [Y = 0.28, Z = 0.02] (dotted 
lines). The first part of the curves shows a well-defined, 
and almost linear, relation between the luminosity and ef- 
fective temperature, with a typical slope d log T e s/d log L 
around -0.13 for [Y = 0.28, Z = 0.02] and -0.11 for 
[Y = 0.25, Z = 0.008], and a very small dependence on 
the initial stellar mass at given metallicity. This feature 
holds as long as the total stellar mass is not significantly 
reduced by stellar winds (i.e. up to the luminosity maxi- 
mum). The subsequent part of the track - corresponding 
to the decline of the luminosity and quick decrease of the 
effective temperature - is controlled by the interplay be- 
tween the weakening of envelope burning and the onset of 
the superwind phase. 

The filled circles show the most luminous AGB stars of 
a sample of long-period variables in the LMC observed by 
Wood et al. (1983) in the near infrared pass-bands JHK. 
The bolometric magnitudes of the observed stars are re- 



Fig. 7. The same as in Fig.6, but with [Y = 0.25, Z = 0.008]. 

scaled to the distance modulus to the LMC, (to — M) = 
18.5, instead of 18.6 used by the authors. The conversion 
of (J — K) colours to effective temperatures is obtained 
with the aid of the relation T cff = 7070/ [(J - K) + 0.88] 
from Bessel et al. (1983). 

Our tracks for [Y — 0.25, Z = 0.008] are consistent 
with the observational data extending in luminosity up to 
about Mboi ~ —7. It is worth recalling that all the stars 
displayed in Fig. 8 are oxygen-rich, in agreement with the 
expectation from envelope burning, preventing the forma- 
tion of carbon stars at high luminosities (Sect. 5). 

4-4- The Mboi — P diagram 

Figure 9 compares in the Afboi — P diagram theoreti- 
cal results for TP-AGB stars with different initial mass 
and composition [Y — 0.25, Z = 0.008] with the obser- 
vational data for long-period AGB variables (Mira and 
OH/IR objects) in the LMC, collected from various au- 
thors. Most of the brightest stars are confined below the 
classical AGB limit, i.e. Mbd ~ —7.1 corresponding to 
the Chandrasckhar mass M c = 1.4M Q as predicted by 
the M c — L relation of Paczynski (1970). Furthermore, 
no theoretical tracks goes beyond this limit. The high- 
est luminosity (Mbd ~ —6.85) is reached by the 5.0M Q , 
[Y = 0.25, Z = 0.008] model, with a M c - 0.95M©. The 
filled triangles indicate some long-period variables which 
are found by Smith et al. (1995) to possess strong Li 
features in their spectra. This is usually interpreted as 
the signature of envelope burning in massive AGB stars 
(Sect. 5.1). The theoretical results seem to be consistent 
with the fact that most of the observed stars more lumi- 
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Fig. 8. Evolutionary tracks in the H-R diagram for stars with 
chemical composition [Y = 0.25, Z = 0.008] (solid lines) and 
[Y = 0.28, Z = 0.02] (dotted lines). The numbers indicate the 
corresponding initial mass in solar units. The filled circles refer 
to the most luminous oxygen-rich AGB stars in the sample of 
long-period variables observed by Wood et al. (1983). 

nous than Mboi ~ — 6 are oxygen-rich. Finally, it is worth 
noticing that the rapid excursion of the tracks towards 
very high periods (P > 800 days) occurs when mass loss 
by stellar winds drastically reduces the envelope mass. In 
the most massive stars (i.e. 4.0, 4.5, 5.OM models), this 
stage coincides with the extinction of envelope burning 
and the recovering of the M c — L relation. 

5. Nucleosynthesis 

5.1. Observational hints for envelope burning 

Envelope burning can heavily alter the surface chemical 
composition of massive TP- AGB stars via the nuclear re- 
actions of the CNO-cycle (and p-p chain), occurring in 
convective conditions. Therefore, the chemical abundances 
observed at the surface of real stars can be used as probes 
of current theories on the structure and evolution of AGB 
models. 

Over the years a considerable amount of spectroscopic 
data on the chemical composition of evolved stars has be- 
come available, thus setting precise constraints on the nu- 
cleosynthesis and mixing processes occurring during the 
evolution of low- and intermediate-mass stars (Smith & 
Lambert 1985; Harris et al. 1985, 1987, 1988; Smith & 
Lambert 1990; Smith et al. 1995). 

Depending on the surface abundance ratio C/O, three 
broad classes of stars can be distinguished: M-giants 



Fig. 9. Bolometric magnitude of long-period variable AGB 
stars as a function of the pulsational period. The data are taken 
from Smith & Lambert 1990, Smith et al. 1995 (triangles); 
Wood et al. 1983, Wood et al. 1992 (squares); and Reid et al. 
1995 (circles). The theoretical tracks correspond to the evolu- 
tion of TP- AGB stars with composition \Y = 0.25, Z = 0.008] 
for a few representative values of the initial mass (indicated in 
Mq by numbers nearby to the corresponding curve). The dot- 
ted lines refer to the results for low-mass stars (M < 3M ), 
already presented in Marigo et al. (1996). 



(C/0< 0.8), S-giants (0.8 < C/0< 1.0), and C-giants 
(C/0> 1.0). According to the present understanding, 
the possible evolutionary sequence along the three classes 
would be the result of the dredge-up of 12 C at successive 
thermal pulses. However, the careful analysis of this sub- 
ject suggests that the real situation is more intrigued and 
interwoven. 

First, a TP- AGB star experiencing envelope burning 
may not become a carbon star (Boothroyd et al. 1993) be- 
cause of the rapid conversion of newly dredged-up 12 C into 
13 C and then into 14 N via the first reactions of the CNO 
cycle. This could offer a viable explanation to the observed 
scarcity of very luminous carbon stars (Mboi < —6.5) in 
the fields of the Galaxy and the fields and clusters of the 
Magellanic Clouds (Westerlund et al. 1991; Costa & Fro- 
gel 1996). Recently, van Loon et al. (1997) have found 
obscured carbon stars in the Magellanic Clouds at lumi- 
nosities up to Mboi ~ —6.8. 

Second, stars with envelope burning are expected to 
rapidly reach high luminosities (Sect. 4.1). This fact would 
speed up the ejection of the envelope by triggering en- 
hanced mass loss and explain the observed paucity of lumi- 
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nous oxygen-rich (M-type) stars in the range —6 > Mboi > 
-7 (Hughes & Wood 1990; Reid et al. 1990). 

However, as far as the analogous situation in clusters of 
the LMC is concerned, there is another effect to be consid- 
ered, i.e. the effect of the past history of cluster formation 
in the age range at which AGB stars are expected to oc- 
cur. In brief, Marigo et al. (1997) have pointed out that 
the lack of these luminous AGB stars in the youngest LMC 
clusters could be ascribed to an epoch of modest cluster 
formation, between ages of 2 x 10 8 and 6 x 10 8 yr ago. 

Third, the detailed information on oxygen and car- 
bon isotopic ratios for galactic giant stars have disclosed 
a more complex scenario (Smith & Lambert 1990) than 
suggested by the classification based on the C/O ratio. C- 
stars can be grouped in two sub-classes according to their 
12 C/ 13 C ratio: 13 C-rich J-type ( 12 C/ 13 C< 10) and N-type 
( 12 C/ 13 C> 10). The low 12 C/ 13 C ratio in J-type stars, 
which is very close to the equilibrium value for the CNO 
cycle, would suggest the occurrence of envelope burning 
(Renzini & Voli 1981). 

In addition to this, extensive spectroscopic searches 
have been carried out in the Magellanic Clouds to detect 
lithium in red giant stars (Smith et al. 1995). The Li I 
6707 A resonance line has been measured in AGB stars, 
mostly belonging to the S-class, within a narrow luminos- 
ity range (—7 Mboi & —6). This feature is commonly 
interpreted as the signature of the Cameron & Fowler 
mechanism for Li production (Cameron & Fowler 1971; 
Sackmann & Boothroyd 1992) taking place in the typical 
envelope burning environment. In brief, 7 Li is produced 
by electron captures on 7 Be nuclei convected from the hot 
regions of the envelope into cooler layers (T < 3 x 1Q 6 K) 
before the reaction 7 Be(p, 7) 8 B proceeds. 

Finally, the expected enrichment in helium and nitro- 
gen caused by envelope burning and mixing episodes is 
usually invoked to interpret the high He/H and N/O ra- 
tios characterizing Planetary Nebulae of Type I (Peimbert 
& Torres-Peimbert 1983). 

5. 2. The nuclear network: our prescriptions 

In this study nuclear burning at the base of the convective 
envelope is followed adopting the instantaneous mixing 
approximation. We check that the lifetimes of the CNO 
nuclei against proton captures are much longer than the 
time scale for convective turnover in the envelope. This 
means that the entire envelope is efficiently mixed through 
the high temperature region at its bottom, so that the new 
nuclei are redistributed all over the envelope faster than 
the time scale over which a significant inhomogeneity can 
be built up (Clayton 1983). At each time, the homoge- 
neous distribution of the nuclei reflects the degree of the 
CNO cycling attained within the burning zones. To give 
an example, a surface 12 C/ 13 C ratio equal or very close 
to 3.4 implies that the entire envelope has been processed 
in the deepest regions at sufficiently high temperatures so 



that the CN cycle (and possibly the complete CNO cycle) 
has reached the state of stationary equilibrium. 

In the case of long-lived nuclei, the reaction rates can 
be properly mass-averaged over the whole convective re- 
gion (Scalo et al. 1975). In contrast in the case of light 
elements such as 7 Be and 7 Li the nuclear lifetimes can be 
comparable to the convective time scale of the envelope. 
In such a case the mass-averaging procedure is not mean- 
ingful and a time-dependent diffusive algorithm should be 
employed in order to properly couple nucleosynthesis and 
mixing. 

In the present study the nuclear reactions rates are 
taken from Caughlan & Fowler (1988). The parameter / 
(0-^1) in their analytical expressions for several uncertain 
rates such as ( 15 N(p, a) 12 C, 17 0(p,a) 14 N, 17 0(p, 7 ) 18 F) 
is set equal to 0.1, as suggested by the authors. In princi- 
ple, the parameter / should be calibrated against chemical 
abundances in real stars. 

Finally, the screening factors are from Graboske et al. 
(1973) and the /3 + -decays are assumed to take place in- 
stantaneously. 

5.3. The third dredge-up 

The analytical treatment of the third dredge-up is de- 
scribed in Marigo et al. (1996), to whom the reader should 
refer for more details. Suffice it to remind that we describe 
this process with the aid of two basic parameters, namely: 

★ M™ ln , the minimum core mass for the third dredge-up 
to occur; 

★ A = AA#drodgo/AM c , the efficiency of the third dredge- 
up, defined as the fraction of the core mass increment, 
AM C , during the previous inter-pulse period, that is 
brought up to the surface (AMd re dge is the dredged- 
up mass). 

Both parameters (M c min = 0.58M©; A = 0.65), have been 
calibrated on the observed luminosity function of carbon 
stars of the LMC (Richer et al. 1979, Blanco et al. 1980, 
Cohen et al. 1981; see also Groenewegen & de Jong 1993 
for similar theoretical results). 

For the sake of simplicity, we do not consider the 
possible dependence of M™ m and A on the stellar mass 
and chemical composition. However, full calculations of 
the third dredge-up actually suggest that this process is 
favoured at lower metallicities and larger envelope masses 
(Wood 1981; Boothroyd & Sackmann 1988c). 

In our model the third dredge-up is considered during 
the TP- AGB phase if two conditions are simultaneously 
met: i) the star has a core with a mass M c > M™ m and, 
ii) the corresponding thermal pulse has already reached 
the full amplitude regime. Moreover, we assume that the 
dredge-up events occur till the end of the TP-AGB evo- 
lution (case A). The adopted chemical composition (in 
mass fraction) of the inter-shell material just before the 
inward penetration of the convective envelope consists of 
4 Hc (76 %), 12 C (22 %) and 16 (2 %) (Boothroyd & 
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Sackmann 1988c). Recently, Herwig ct al. (1997) have 
shown that the introduction of additional partial mixing 
in the treatment of the third dredge- up may lead to signif- 
icantly changed inter-shell abundances (e.g. 4 He (25 %), 
12 C (50 %) and 16 (25%)). 

Our simple treatment of the third dredge-up is mainly 
justified by the unsatisfactory understanding of the third 
dredgc-up phenomenon even in stellar models that care- 
fully follow thermal pulses. On one hand most authors 
still do not find significant dredge-up in their calculations 
(Vassiliadis & Wood 1993; Wagcnhuber 1996), or assume 
the efficiency A as an input quantity to be fixed by suitable 
calibrations (Forestini & Charbonnel 1997). On the other 
hand, the difficulty of too poor a dredge-up in low-mass 
stars seem to be overcome to some extent (Boothroyd & 
Sackmann 1988c; Frost & Lattanzio 1996; Straniero et al. 
1997, Herwig et al. 1997), but extensive calculations over 
a wider range of stellar masses and metallicities are re- 
quired in order to derive quantitative predictions which 
can be readily used in analytical models. 




Time (1 3 vr) 



Fig. 10. Evolution of the surface abundances (by mass) of C, 
N, O and their ratios for the 5.0M Q , [Y = 0.25, Z = 0.008] 
TP-AGB star. The solid lines refer to case A, whereas the 
dashed lines describe the final freezing of the abundances in 
case B. See the text for further explanation. 

In Fig. 10 we plot the time evolution of the surface 
abundances of C, N, and O for the 5.OM , [Y = 0.25, Z = 
0.008] star, starting from the onset of the thermally puls- 
ing regime till the complete ejection of the envelope. The 
saw-like trend of the C abundance is caused by the com- 
bined effect of the dredge-up episodes at each flash - 
clearly visible in the sudden spikes -, and envelope burn- 



ing during the subsequent inter-flash period, which rapidly 
converts C into N. 

Note that for most of the lifetime a TP-AGB star is 
expected be to oxygen-rich, in agreement with the obser- 
vational indication of the most luminous AGB stars (see 
Sects. 4.4 and 5.1). In this case, the transition to a carbon 
star eventually happens when envelope burning is extin- 
guished, and 12 C enrichment due to convective dredge-up 
becomes dominant. Very recent observations (van Loon et 
al. 1997) in N-band and near- infrared photometry, aimed 
to detect dust-enshrouded AGB stars in the LMC, indi- 
cate that very luminous carbon stars (Mboi ~ —6.8) are 
quite rare, but not absent. 

Actually, spectroscopic estimates of the abundances in 
Planetary Nebulae (PN) indicate that high values of the 
ratios (N/O > 0.5) and (He/H > 0.125) typical of Type I 
PN are associated with rather low C/O ratios (in general 
C/O < 1). 

To investigate the effect of a different prescription we 
draw in Fig. 10 the dashed horizontal lines corresponding 
to the situation of a concomitant shut-down of the third 
dredge- up when envelope burning extinguishes (case B). 
In this case the formation of the carbon star is inhibited, 
and the final value of the C/O ratio is approximately 0.55, 
closer to the value in Type I PN. 

These simple experiments are interesting in view of 
exploring different effects, without leading to draw defini- 
tive conclusions. Indeed, the competition between dredge- 
up and envelope burning is a crucial point that has not 
yet been solved by full AGB calculations, due to the 
uncertainties in the treatment of convection (as already 
mentioned), and additionally on the envelope mass, i.e. 
on mass-loss. Therefore, the critical envelope masses for 
the shut-down of dredge-up and envelope burning, re- 
spectively, are not known exactly up to now. In the case 
of the 6Mq,Z = 0.02 model of Lattanzio & Boothroyd 
(1997) envelope burning seems to extinguish earlier than 
the dredge-up episodes do. Our case A resembles this sit- 
uation. 

For each model star in Table 1 we report the critical en- 
velope mass, M e ™*, at which hot-bottom burning is found 
to be no longer operative (e.g. when the contribution to 
the total luminosity falls below 1%) and the number of the 
dredge-up episodes occurred beyond this point according 
to case A. 

5.4- Type I Planetary Nebulae: how to get high He/H ? 

Chemical abundances of PN in the Galaxy and Magellanic 
Clouds are an important observational constraint on pre- 
scriptions for nucleosynthesis, mixing, and mass-loss used 
to calculate the evolution of low- and intermediate-mass 
stars. 

In the following, we will limit ourselves to address the 
question of the high overabundances of helium (He/H> 
0.125) and nitrogen (N/0> 0.5) detected in Type I PN. 



Envelope Burning 



13 



Table 1. Envelope mass, at the extinction of envelope 

burning (when present) and corresponding inter-pulse number, 
Np. This stage is defined as soon as Leb/L < 0.01. ANp in- 
dicates how many dredge-up events happen beyond this point 
till the end of calculations (case A). 



z 


M(M ) 


M^(M Q ) 


N P 


AiVp 


0.008 


4.0 


2.417 


117 


57 




4.5 


2.365 


86 


60 




5.0 


2.221 


59 


70 



data. It is worth noticing that a successful situation might 
occur if only few dredge-up events (contrary to our case 
A, see Table 1) follow the cessation of envelope burning, 
so that C/O could remain smaller than unity. 

(3) Under the assumptions of case B (i.e. the dredge- 
up episodes cease when envelope burning extinguishes), 
we investigate the effect of increasing the efficiency of the 
third dredge-up (e.g. setting A = 0.9). This actually leads 
to higher helium abundances, but always well below the 
observational value of He/H> 0.15. 



Z M(Mq) Ai£|*(M ) Np ANp 

0.02 4.0 no envelope burning 

4.5 3.053 15 41 

5.0 2.679 30 53 



According to the present understanding (Renzini & Voli 
1981; Groenewegen & de Jong 1994), high values of ni- 
trogen would mainly be the result of envelope burning in 
massive AGB stars (the contribution from the first and 
second dredge- up being marginal). 

This indication is confirmed by our results (case B) 
shown in the top panel of Fig. 11, where the N/O ratio 
attained after the second dredge-up eventually increases 
to considerable values during the TP-AGB evolution in 
partial agreement with the observations. In the case of the 
4.5, and 5.0M©, [Y = 0.25, Z = 0.008] stars the predicted 
N/O ratios exceed, roughly by a factor of two, the upper 
boundary (log N/O ~ 0.5) set by measurements of chem- 
ical abundances in PN of the LMC. The corresponding 
final C/O ratios (bottom-panel of Fig. 11) are also consis- 
tent with observations. The low values (C/0< 1) are the 
result of envelope burning rapidly converts carbon into 
nitrogen. In addition to this, the observational data seem 
also to suggest that during the final stages of the TP-AGB 
evolution the third dredge-up is likely to vanish, thus keep- 
ing low the C/O ratio (see Fig. (10)). 

As far as the helium content is concerned, one can im- 
mediately notice in Fig. 11 that none of our results extends 
beyond He/H ~ 0.14, whereas the observational estimates 
can be as high as ~ 0.2. 

Several experiments are performed to understand the 
reason of the marginal disagreement: 

(1) With the present prescriptions a more efficient en- 
velope burning (for instance obtained by increasing the 
mixing length parameter) can be ruled out as a larger 
production of helium also implies a significant increase of 
nitrogen, thus easily exceeding the upper limit on the N/O 
ratio indicated by the observations. 

(2) High values of He/H could be obtained by letting 
the third dredge-up proceed till the very end of the TP- 
AGB phase (case A). In such a case, however, ratios C/O 
is expected to get » 1, in contrast with the observational 
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Fig. 11. N/O ratio (top-panel) and C/O ratio (bot- 
tom-panel) as a function of He/H for the 4.0, 4.5, 5.0M©, 
[Y = 0.25, Z — 0.008] stars calculated according to the pre- 
scription of case B. The abundances are expressed as fractions 
by number. The three lines connecting the series of symbols 
(one for each value of the stellar mass as indicated) correspond 
to three theoretical abundances: 1) after the first dredge-up, 
2) after the second dredge-up, 3) mass-averaged in the ejecta 
during the last 5 x 10 4 yr before the end of the AGB phase. 
The asterisks refer to the initial abundance ratios. The dashed 
areas roughly correspond to the observed domain of type-I PNe 
in the LMC. 

(4) Finally, we explore the possibility of a much deeper 
penetration of the convective envelope during the sec- 
ond dredge-up. Let us consider the case of the 5.0M©, 
[Y = 0.25, Z = 0.008] star. The second dredge-up makes 
the surface ratio (by number) He/H to increase from 
<~ 0.092 up to ~ 0.124. Correspondingly, the hydrogen 
Yh an d helium Yho abundances (by number) change from 
(Y& ~ 0.726, Y^ c ~ 0.066) to (Y£ ~ 0.662, Y^ c ~ 0.082), 
where the superscripts 1 and 2 refer to the values just 
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before and after the second dredge-up, respectively. The 
dredged-up mass (in M Q ) is given by: 



AM 2D = M(QxY - Q^nv) 



(19) 



where M (in M Q ) is the total stellar mass, Qxy ~ 
0.252853 defines the X — Y discontinuity when the H- 
burning shell extinguishes, and Q™„v ~ 0.186333 corre- 
sponds to the maximum penetration of the envelope con- 
vection during the dredge-up, with Q — M r /M (mass co- 
ordinate). In this case, AM 2D amounts to 0.3326M Q . 

With aid of the above quantities, we can evaluate 
the maximum amount of material, AM 2 ™ that can be 
dredged-up if the envelope extends down to the upper 
boundary (Qhc ~ 0.184678) of the He-burning shell, 
which is an insuperable barrier to the inward displacement 
of the convection. The dredged-up material is 



AM 2 7 = M(Qxy ~ Qhc) 



(20) 



which gives AM™ = O.341OM . Note that this value 
slightly exceeds the actual value already calculated with 
Eq. (19). 

For the sake of comparison, we estimate also the 
amount of mass, AM' 2 d, that should be dredged-up in or- 
der to get the typically high values of the helium content 
measured in Type-I PN. As an example, let us consider 
the case of He/H~ 0.176, or equivalently Y$ ~ 0.582, 
Y* e ~ 0.1025. Denoting by Y$ e ~ (1 - Z)/A = 0.248 the 
helium abundance (by number) in the dredge-up mate- 
rial, and M env = M(l — Qxy) ~ 3.73M Q the mass of the 
convective envelope just before the penetration, we have: 



yi 

AM',n = Mr. Hc 



Y 2 ' 



Y 



He 



yd 



(21) 



Equation (21) yields AM' 2D ~ 0.93M©, which is 
higher than the actual value by a factor of almost three 
and, above all, it largely exceeds the maximum value al- 
lowed, AM 2 ™ ax . 

Similar estimates, carried out for the other star masses 
considered in this paper, lead us to the conclusion that in 
order to get high values of the He/H ratio (> 0.16) the 
second dredge-up cannot play the prevalent role. 

All things considered, we expect that the dominant 
causes of helium enrichment are nucleosynthesis and mix- 
ing occurring during the TP-AGB phase and, mainly, cu- 
mulative effects of recurrent dredge-up episodes. 

5. 5. Secondary and primary synthesis of CNO nuclei 

During the TP-AGB evolution of stars undergoing enve- 
lope burning, at each dredge-up event some primary 12 C 
and 16 O nuclei (produced by a-capture reactions) are in- 
jected into the envelope where they can partake in the 
CNO cycle for the conversion of hydrogen into helium. 
Therefore, as a consequence of chemical pollution by the 



third dredge-up and subsequent nuclear burning at the 
base of the envelope, the surface distribution of the CNO 
nuclei consists of two components: 

— A secondary component produced by nuclear reactions 
involving only seed metals originally present in the star 
at the epoch of its formation; 

— A primary component synthesized by a chain of nu- 
clear burnings beginning from H and He and hence 
independent of the original metal content. 

The surface abundance (by number) Yi of any element 
i of the CNO group at time t, can be expressed as Yi = 
Y i s +Y i p , where Yf and Y? are the secondary and primary 
fractions, respectively. At the time t = t + dt, because of 
nuclear burning at the base of the convective envelope, 
the element i will achieve a new surface abundance Y i = 
Y? + Y? , with a different relative partition of the two 
components. 

The method to follow the time evolution of both the 
secondary and primary nuclear synthesis is quite simple, 
in virtue of the following considerations. The generic form 
of the differential equation describing the time variation 
of the abundance Yi can be written as: 



dt 



(22) 



or equivalently, 
^L+Xll = ^ (± )^.y H( y 4 S + yP) (23 ) 



where the summation is performed over all the nuclear 
reactions that create the element (sign +) and destroy 
it (sign — ), and the various /3-decays are neglected. Two 
important points are to be noticed: i) the nuclear rates 
Rij are functions of temperature and density only, and 
ii) the CNO nuclei react only with protons and not 
among each other. This latter condition is clearly ex- 
pressed by the product YgYi in Eq. (22), where y H is 
the hydrogen abundance. Since the only possible reac- 
tions are proton-captures on CNO nuclei it follows that in 
Eq. (23) mixed products of secondary and primary abun- 
dances (e.g. Yp"Y?) are not possible. Therefore, Eq. (23) 
can be properly split into two independent equations: 



dY 



3 



dt 
dY? 



(24) 
(25) 



which can be integrated separately over the time dt to 
give the new secondary Y t s and primary Y? fractions, 
such that Y? + Y? = Y i . In other words, the differential 
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Table 2. Chemical abundances (in mass fraction) in stars with initial mass Mi and original metallicity Z = 0.008 at three 
stages, namely: i) at the zero-age main sequence (Initial); ii) soon after the first dredge-up (1 D); and iii) soon after the second 
dredge- up (2D). 







H 


3 He 


4 He 


12 C 


i3 C 


14 N 


15 N 


16 


17 o 


18Q 


Initial 


0.7420 


2.675E-05 


0.2500 


1.371E-03 


1.652E-05 


4.238E-04 


1.672E-06 


3.851E-03 


1.560E-06 


8.680E-06 


1 D 


Mi 
























4.0 


0.7417 


1.125E-04 


0.2503 


1.167E-03 


5.927E-05 


6.167E-04 


1.055E-06 


3.850E-03 


1.607E-06 


8.331E-06 




4.5 


0.7417 


9.285E-05 


0.2502 


1.160E-03 


5.729E-05 


6.263E-04 


1.044E-06 


3.850E-03 


1.631E-06 


8.249E-06 




5.0 


0.7259 


6.845E-05 


0.2661 


9.833E-04 


4.765E-05 


1.098E-03 


9.139E-07 


3.559E-03 


3.196E-06 


7.003E-06 


2 D 


Mi 
























4.0 


0.7032 


9.179E-05 


0.2887 


9.164E-04 


4.959E-05 


1.261E-03 


8.402E-07 


3.459E-03 


4.623E-06 


6.595E-06 




4.5 


0.6809 


7.394E-05 


0.3111 


9.039E-04 


4.871E-05 


1.379E-03 


8.238E-07 


3.342E-03 


4.197E-06 


6.445E-06 




5.0 


0.6625 


5.991E-05 


0.3295 


8.756E-04 


4.654E-05 


1.502E-03 


8.052E-07 


3.243E-03 


3.691E-06 


6.219E-06 



Table 3. The same as in Table 2, but with initial metallicity Z = 0.02. 







H 


3 He 


4 He 


i2 C 


13 C 


14 N 


15 N 


16q 


17 o 


18 


Initial 


0.7000 


2.996E-05 


0.2800 


3.427E-03 


4.130E-05 


1.060E-03 


4.180E-06 


9.626E-03 


3.900E-06 


2.170E-05 


1 D 


Mi 
























4.0 


0.6778 


8.915E-05 


0.3021 


2.420E-03 


1.222E-04 


2.687E-03 


2.262E-06 


9.003E-03 


1.599E-05 


1.749E-05 




4.5 


0.6776 


7.410E-05 


0.3022 


2.444E-03 


1.254E-04 


2.690E-03 


2.253E-06 


8.966E-03 


1.358E-05 


1.764E-05 




5.0 


0.6986 


6.937E-05 


0.2814 


2.692E-03 


1.372E-04 


1.885E-03 


2.408E-06 


9.546E-03 


8.416E-06 


1.937E-05 


2 D 


M 
























4.0 


0.6713 


8.684E-05 


0.3087 


2.363E-03 


1.221E-04 


2.839E-03 


2.197E-06 


8.906E-03 


1.738E-05 


1.713E-05 




4.5 


0.6485 


6.961E-05 


0.3314 


2.305E-03 


1.230E-04 


3.189E-03 


2.123E-06 


8.586E-03 


1.446E-05 


1.665E-05 




5.0 


0.6307 


5.830E-05 


0.3494 


2.282E-03 


1.241E-04 


3.432E-03 


2.072E-06 


8.341E-03 


1.218E-05 


1.633E-05 



equation (23) can be seen as a linear transformation 1 of 
the variables Yj S and Yf '. 

In Figs. 12 and 13 the evolution of secondary and pri- 
mary 12 C, 13 C, 14 N, 15 N, 16 0, 17 0, 18 surface abun- 
dances (by mass) is plotted as a function of time for the 
4.5M TP-AGB star with [Y = 0.25, Z = 0.008]. It can 
be noticed that for a few CNO nuclei (e.g. 12 C, 13 C, 14 N, 
and 15 N), the primary component (which is initially equal 
to zero by definition) of the surface abundance grows until 
it overcomes the secondary one. 



ical compositions are 4.0, 4.5, 5.0M Q and [Y — 0.28, Z = 
0.02] and [Y = 0.25, Z = 0.008]. These data complete 
the sets of chemical yields from lower mass stars (0.7 < 
M/Mq < 4.0) already presented by Marigo et al. (1996). 
The "old" results for the 4.0M© models in Marigo et 
al. (1996) should be replaced with the "new" results of 
this study. For each elemental species i we calculate the 
amount of mass M y {i) which is newly synthesized by nu- 
clear burning in the star and returned to the interstellar 
medium during the entire lifetime (T): 



6. Yields of chemical elements 

Tables 4 and 5 give the chemical yields from stars which 
experience envelope burning. The initial masses and chem- 

1 We recall that a function f(x, y) is a linear transformation 
if f{ax + by) = af{x) + bf(y). 



M y (i) = j\ X {i)-X\i)]^&t (26) 

where X(i) and X (i) are the current and the initial sur- 
face abundance, respectively; dM/dt is the current mass- 
loss rate. 



16 



Envelope Burning 




2 3 
Time (10 6 yr) 



Fig. 12. Secondary and primary surface abundances (by mass) 
of 12 C, 13 C, and 14 N as a function of time for the 4.5M©, 
[Y = 0.25, Z = 0.008] star in the TP-AGB phase (case B). 
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Fig. 13. The same as in Fig. 12, but for the surface abundances 
of 15 N, 16 0, 17 Q, and 18 (case B). 



It is worth recalling that with the adopted mass-loss 
rates along the TP-AGB phase (Vassiliadis & Wood 1993), 
the reduction of stellar mass during the E-AGB is indeed 
negligible. 

As far as the CNO nuclei are concerned, the entry 
M y (CNO) in Table 4 corresponds to the total yield of all 
the isotopes, whereas the entries of Table 5 explicitly gives 
both the secondary and primary yield for each element. 
These are calculated with: 



MS(i) = J o T [X s (i)-X\i)]^dt 
f T vT >, ,dM , 



My (i) 



(27) 



(28) 



The total yield of any element is simply M y = M y (i) + 
M y p ( ? ). 

According to the above definitions, the primary yields 
can be only > 0, whereas secondary yields can be either 
> or < 0, in the respective cases that the mass-averaged 
abundance of the element in the ejecta is greater, equal 
or smaller than its original value at the formation of the 
star. 

The reader should refer also to Tables 2 and 3 for an 
easier understanding of the effects on the surface abun- 
dances caused by different mixing episodes during the evo- 
lution of the star. Suffice it to recall that the both the first 
and second dredge-up affect (by increasing or decreasing) 
only the secondary components of the CNO surface abun- 
dances. In fact, in these cases the envelope is polluted by 
material which has undergone CNO-cycling of nuclei syn- 
thesized from metal seeds originally present in the star. 

From a general analysis of the total yields and their 
components (primary versus secondary) it turns out that 
three alternatives are possible: 

1. M y s (i) > and M y p (i) > so that M y > 

2. AfS(i) < and M y p (i) > so that M y > 

3. M y s (i) < and M*(i) > so that M y < 

The first case applies, for instance, to the 13 C abun- 
dance which is significantly enhanced by the first dredge- 
up (secondary origin), and subsequently by envelope burn- 
ing, thanks to proton-capture reactions on 12 C nuclei (sec- 
ondary and primary nucleosynthesis). It results that the 
primary yield of 13 C exceeds (by few factors, depending 
on the efficiency of envelope burning) the secondary con- 
tribution in all the stars under consideration, but for the 
4.0M Q , [Y = 0.28, Z = 0.02] star, in which the effect of the 
first dredge-up prevails. Likewise, the relative weight of 
the secondary and primary components on the total yield 
of 14 N depends on the interplay between the secondary en- 
richment due to the first and second dredge-up episodes, 
and the primary contribute given by the CNO-cycle during 
envelope burning. The highest value of My( 14 N) is pro- 
duced by the 4.0M Q , [Y = 0.25, Z = 0.008] star, in which 
envelope burning has been active for the longest time. 
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Table 4. Total yields (in M ) from intermediate-mass stars with initial metallicity Z, and mass M\ (case A). The amount of 
mass, M e j, ejected over the entire evolution, and the final mass of the core, Mf, are indicated in solar units. 



z 


Mi 




Mf 


M y (H) 


71 ir /3 TT \ 

M y { He) 


Tt jt /4tt \ 

M y { He) 


My (UNO) 


0.008 


4.0 


3.059 


0.941 


-2.218E-01 


-8.130E-05 


1.877E-01 


3.386E-02 




4.5 


3.545 


0.955 


-3.374E-01 


-9.478E-05 


3.138E-01 


2.363E-02 




5.0 


4.018 


0.982 


-4.024E-01 


-1.075E-04 


3.885E-01 


1.453E-02 


0.02 


4.0 


3.143 


0.857 


-1.121E-01 


1.594E-04 


1.087E-01 


8.798E-03 




4.5 


3.604 


0.896 


-2.178E-01 


-7.890E-05 


2.059E-01 


1.202E-02 




5.0 


4.061 


0.939 


-3.257E-01 


-1.188E-04 


3.129E-01 


1.379E-02 



Table 5. Total yields from intermediate-mass stars (case A). The notation is the same as in Table 4. For each isotope of the 
CNO group, the secondary (S) and primary (P) yields are distinguished. 



z 


M 


M cj 


M f 


S/P 


M y ( 12 C) 


M y ( 13 C) 


M y ( 14 N) 


M y ( 15 N) 


M y ( ie O) 


M y ( 17 0) 


M y ( ls O) 


0.008 


4.0 


3.059 


0.941 


S 


-3.884E-03 


3.546E-05 


5.540E-03 


-4.905E-06 


-2.735E-03 


8.199E-04 


-2.645E-05 










P 


8.014E-03 


4.845E-04 


2.307E-02 


7.100E-07 


2.422E-03 


1.246E-04 


1.532E-08 




4.5 


3.545 


0.955 


s 


-4.336E-03 


9.594E-05 


1.195E-02 


-5.532E-06 


-9.652E-03 


1.179E-03 


-3.068E-05 










p 


7.741E-03 


3.601E-04 


1.475E-02 


4.369E-07 


1.452E-03 


1.333E-04 


1.157E-08 




5.0 


4.018 


0.982 


s 


-4.855E-03 


1.266E-04 


1.377E-02 


-6.281E-06 


-1.099E-02 


1.299E-03 


-3.481E-05 










p 


7.177E-03 


2.439E-04 


6.697E-03 


1.947E-07 


1.046E-03 


5.667E-05 


3.951E-09 


0.02 


4.0 


3.143 


0.857 


s 


-3.503E-03 


3.298E-04 


5.516E-03 


-1.253E-05 


-2.556E-03 


4.230E-05 


-1.834E-05 










p 


8.235E-03 


1.411E-05 


3.327E-07 


8.870E-12 


7.498E-04 


1.689E-09 


1.797E-14 




4.5 


3.604 


0.896 


s 


-1.090E-02 


2.801E-04 


1.545E-02 


-1.447E-05 


-4.414E-03 


3.098E-04 


-7.802E-05 










p 


8.743E-03 


3.028E-04 


1.415E-03 


4.364E-08 


9.299E-04 


9.966E-07 


1.068E-10 




5.0 


4.061 


0.939 


s 


-1.263E-02 


1.879E-04 


1.980E-02 


-1.628E-05 


-9.853E-03 


3.293E-03 


-8.791E-05 










p 


8.345E-03 


2.184E-04 


3.489E-03 


1.035E-07 


1.033E-03 


1.469E-05 


9.511E-10 



It is worth remarking that intermediate-mass stars with 
M > A.QMq are a source of primary 14 N for the chemical 
enrichment of the interstellar medium. These stars could 
play an important role in view of interpreting, with the 
aid of chemical evolutionary models of galaxies, the ob- 
served trend of the data in the log(N/0) vs. log(0/H) 
plane, which suggests the existence of a significant pri- 
mary component of 14 N (see Matteucci 1997). A positive 
yield is predicted also for 17 0, thanks to the first dredge- 
up which increases the secondary surface abundance, and 
to CNO-cycling of primary nuclei during envelope burn- 
ing. 

The second case applies to 12 C, for which the surface 
abundance is first decreased as a consequence of both the 
first and second dredgc-up episodes, and then increased 
thanks to the convective dredge-up at thermal pulses. In 
all the stars considered, the positive net yield M y ( 12 C) is 
of primary origin only. 



Finally, the third case refers to the other CNO iso- 
topes ( 15 N, 16 0, and 17 0), for which the negative contri- 
bution of the secondary component exceeds the positive 
one of primary nature. It follows that for these elements 
no chemical enrichment of the interstellar medium is ex- 
pected from the stars under consideration. We notice that 
as far as 16 is concerned, My ( 16 0) crucially depends on 
the number of dredge-up episodes that occurred during 
the TP-AGB phase. The 4.0M Q , [Y = 0.25, Z = 0.008] 
model star - which experiences the greatest number of 
thermal pulses - produces a primary yield nearly as high 
as the secondary one. 

7. Concluding remarks 

In this study we have developed an original method to 
follow the TP-AGB evolution of stars with envelope burn- 
ing. To this aim we formulated a suitable model for static 
envelopes, in which the equation of energy balance is in- 
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eluded. The structure of the envelope is then solved in 
details with the aid of two boundary conditions having a 
general validity, so that the new algorithm can be applied 
also to the case with no envelope burning. The method we 
are proposing is very flexible and it can be employed to 
investigate different questions concerning the occurrence 
of envelope burning. 

Among other results, we draw attention on the prob- 
lem raised by the high values of helium content in Typc-I 
PN, which cannot find an easy explanation without violat- 
ing the observational abundances of nitrogen. This is an 
important conclusion setting strong constraints not only 
on nucleosynthesis, but also on other evolutionary aspects 
of AGB stars, e.g. mass-loss, M c — L relation, lifetimes, 
distribution on the H-R diagram, and pulsational proper- 
ties. 
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